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Abstract 

Cross sections of the e + e~ annihilation into one photon plus two vector mesons, such as J /ipp, 
J/ip(f>, and p4> at the center of mass energy y/s = 10.58 GeV are calculated. These shed light on 
the measurement of the e + e~ annihilation cross section via radiative return in i?-factories. At 
another center of mass energy yfs = 3.097 GeV, namely the J/ip resonance peak, the processes 
e + e~ — > jVV and e + e~ — > 77^ (with V being p, lo, or cf)) are calculated with similar method. 
These calculations give estimation of the background levels in the study of radiative or double 
radiative decays of J/tp. 
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I. INTRODUCTION 



The measurement of the e + e~ annihilation cross section via initial state radiation (ISR) 
turns out to be very fruitful at the £>-factories due to the large data samples collected on the 
T(AS) resonance. The study of the final states with a charmonium (J/ift or ip(2S)) and light 
hadrons such as 7r + 7r~ and K + K~ results in the discovery of many unexpected charmonium- 
like states, such as the F(4008), F(4260), F(4360), and y(4660) 0, 0, 0, 3 . But one should 
notice that, e + e~ — > ^ypJ/ip is actually an important background in the measurement of the 
7T + 7T~J/ip production via ISR processes [l|, y[; and this is also true in the case when final 
state 7r + 7T~il)(2S) was concerned @, 0]- Another similar case is in the study of K + K~ J/ift via 
ISR, here background from ^(ftj/ift should be considered. Actually, of a recent observation 
at Belle |5(, this background has been estimated and subtracted by a requirement on the 
K + K~ invariant mass does not agree with a (ft. While the subtraction of ^(ftJ/ift background 
is possible due to the narrow width of the (ft meson, it is almost impossible for the ir + TT~J/ip 
and 7r + 7r~if)(2S) case, where the ir + 7T~ invariant mass distribution is not very different from 
the wide p resonance, so a calculation of the cross section is desired. 

At lower energies, in the studies of J/ift -> ^VV at BES @, 0, 8, where V denotes 
light neutral vector meson such as p (or ou, (ft), direct productions e + e~ —>■ 7VV are back- 
grounds too. Analogously, e + e~ — > •yyV are backgrounds in the J/ift double radiative decay 
channels [s[ as well as in the energy region of few tens of MeV around (ft resonance where 
processes e + e~ — > um° and e + e~ — > pir are measured at KLOE 111!. We can see that, 



although the direct production of e + e~ — > jVV or e + e~ — > 77^ is believed only at order 
0(a 3 ), the high luminosity at the (ft, 73-factories and the forthcoming BESIII will provide 
opportunities to explore these rare processes. 

The process e + e~ — > hadrons at center of mass (CM) energy y/s far below the Z° mass 
is dominated by annihilation via a single virtual photon with charge- conjugation parity 



C = —1. Recently BaBar [12| presented the first observation of the exclusive reactions 
e + e~ — > p°p° and e + e~ (ftp , in which the final states are even under charge conjugate, and 
therefore cannot be produced via just a single photon. A possible interpretation is that these 
C = +1 final states are produced in the two-virtual-photon annihilation processes, i.e., they 
arise from e + e~ annihilation into two virtual photons with each virtual photon converts into a 
vector meson. The rates predicted by the above mechanism can be computed unambiguously 
using the effective vector meson-photon couplings determined from the leptonic widths of the 



mesons. Based on this assumption, the authors of Ref. [13| calculated the cross sections of 
series of these processes and the results are in good agreement with BaBar's measurements. 

In this article, we calculate the cross sections of several e + e~ — > jVV processes, such as 
e + e~ annihilate into one photon plus J/ipp, J/ift<f) an d P<ft at the CM energy y/s = 10.58 GeV 
using the same method proposed in Ref. [13j |. We also calculate e + e~ annihilate into one 
photon plus VV, where VV denotes ojcft or pcft at y/s = 3.097 GeV, in addition to calculate 
e + e~ annihilate into two photons plus a p, a (ft or an u. The structure of this paper is, after 
a short introduction of the model, we give some details of how to calculate the amplitude of 
e + e~ — > (71)7(771)7*, then the numerical results and finally a brief discussion. 



II. AMPLITUDE OF e+e" -> (71)7(771)7* 

Our motivation focuses on investigating the non-resonance contributions to the processes 
of e + e~ — > 7V1V2 and 77^, where the vector mesons of V\ and V2 have different quark 
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contents. We assume the generic reaction e + e~ — > 7V1V2 proceeds via an intermediate 
e + e~ — > 77*7* process, and the two virtual photons converting into V\ and V2 with effective 
couplings e/ fx and ej respectively. Then the differential cross section simply reads as 



e 
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where dcr 77 * 7 * (m^ , m^) is the differential cross section of e + e~ — > 77*7*, which depends on 
the virtual photon masses, i.e. the masses of final vector mesons. The effective photon-meson 
couplings can be directly defined using the leptonic widths of the vectors 

(2) 

when narrow widths approximation is used. Notice that the p meson cannot be described 
as a narrow vector meson properly for its somewhat large width, a complete consideration 
should take an integral over its mass distribution. But from a practical point of view the 
narrow width approximation works well, and the p meson mass distribution only contribute 




a correction less than 8% [131 ] . So we will adopt narrow widths approximation in following 
computation. For another similar reaction e + e~ — > 77^, the cross section formula is similar 
to Eq. ([1]) except one virtual photon be replaced by a real one as long as the interference 
between the two real photons can be neglected. 

In order to calculate the cross sections of e + e~ — > ^VV or e + e~ — > j'jV, from Eq. (TjQ) 
it's clear that an essential work is to calculate the cross section of a pure QED process, 
which can be represented by a more general form e + e~ — > (71)7(771)7*, i.e. an electron and 
a positron annihilate into n real and m virtual photons. Although the amplitudes of these 
processes can be derived from corresponding Feynman diagrams unambiguously, the length 
of the calculation formulae should increase very fast and become very tedious even when 
n + m = 3, let alone the finial states containing more than three photons. As the best to 
our knowledge, present calculations remain on the electron positron annihilation into three 
real photons [15[ or two photons up to the next to leading order [ly] . In practice we utilize 
two specific packages, FeynArts [TtJ and FeynCalc within the symbol calculation tool 
Mathematica to do the programme work to overcome the complexity problem we mentioned 
above and to avoid the potential mistakes which may arise from the lengthy formulae. As 
our interest lies in those particular processes whose final states involving virtual photons, we 
introduce a new model package which contains additional time-like massive photons obey 
the same dynamics of the standard QED photons and should be considered as an extension 
to the built-in QED model. The Feynman diagrams of a relatively simple process e + e~ 
annihilate into one real photon and two virtual photons are depicted in Fig. [1] at tree order 
as an example. 

For a general process e~(pi)e + (p2) — * 7(^1)7(^2)7(^3), the starting point is the scattering 
formula LU 



da(a -+/3) = dT(a -> /3)/$ Q = (2vr) 4 ^ 1 |M^| 2 5 4 (p /3 - Pa )d(3 (3) 

with 



u, 



N(|±^). (4) 
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FIG. 1: Feynman diagrams of e + e — > 77*7* at tree order. 



Here we get u = 2 when the electron and positron masses are neglected, d/3 is the 3-body 
phase space. The corresponding amplitude 



F-M Pa = v{p 2 )f{h 



\k 2 



s (A; 3 )w(pi) + {perm(A;i,fc2,A;3)} (5) 



can be read from FeynArts directly with our modified-QED model, where u, v are spinors 
of initial states, e's are polarization vectors of final states, Pi and fcj are corresponding 
momenta, perm(ki, k 2 , k 3 ) means all the possible permutation of kx, k 2 , and k 3 . And 



F = (27r) 9//2 y 8E[E 2 E 3 is a scalar factor arising from the representation differences between 
FeynArts and FeynCalc packages. Here we use the standard Lorentz covariant formula de- 
rived from FeynArts instead of adopting the helicity form which is proposed in Ref. 13 
just for the technical reason, because in Ref. [2(| polarization vectors, fermion spinors, and 
amplitude matrices are all expressed in explicit forms. In order to speed up the numerical 
integral, we derive a new phase space formula in Lorentz covariant form: 



|fci| 2 |fe 3 |sin^d^id^3d0id^d^ _ |fc 1 | 2 |fc 3 |dft 1 d0 3 d.E^d.E^ 



(6) 



sin 9 2 sin(0 2 - <Ps) sin 9 2 sin(0 2 - 03 ) 

where the sum over all possible solutions of energy-momentum conservation respect to <p 2 
and $3 is implied. The formula is derived with the transformation of 5 function: 



8{h{e u e 2 ))8{j 2 {e u e 2 )) = -61)5(6, -el) a- 1 

where fi and f 2 are general functions depending on 61 and 6 2 , 

9A 



e 1= eo,e 2 =e° 



A = det 



dfr 

de 1 



9A 

dh 

de 2 



(7) 



(8) 
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and 9® and 9% are the solutions of equations fi(9i,9 2 ) = and f2(0i,@2) = 0, means 
a sum over all possible solutions. Considering the energy-momentum is conserved and the 
final state particles are on-shell, we obtain Eq. (jSj). Combining Eq. (j3J) with Eqs. (jlj), (JSj), 
and (jSD, we obtain the cross section formula used in our program [21 . 



III. NUMERICAL RESULTS 

With the amplitudes of e + e~ — ■> 77*7* and e + e~ — > 777*, we obtain the cross sections 
of e + e~ — > 7VV and e + e~ — > 77V by Eq. ([I]). Since the final analytic results are extremely 
lengthy, we only provide the numerical results here for compactness. All the parameters 
used in our computation without explicit exception are quoted from PDG [221 ]. We show 
the cross sections of processes e + e~ — > jVV at y/s = 10.58 GeV in Table HJ with VV being 
J/i[>p, J/il>4> and pcf) when the invariant mass of the meson pair ranging from threshold to 
5.2 GeV/c 2 . As examples for illustration, the differential cross section of e + e~ — > jJ/ipp 
versus the invariant mass of J/ipp is drawn in Fig.[2j Similarly, the cross sections of processes 
e + e _ — ► 7VV at y/s = 3.097 GeV are exhibited in Table HT1 this time VV represents <puj or 
<fip, with the invariant masses ranging from their thresholds to 3.0 GeV/c 2 . Figure [3] shows 
the differential cross section of e + e~ — > 70p versus the invariant mass of <pp as an example. 
Finally, we provide the cross sections of e + e~ — > 77V at y/s = 3.097 GeV in Table IHIl 
here V represents a single meson which is p, <p or with the invariant mass of one photon 
and one meson varying from their thresholds to 3.0 GeV/c 2 . Here only the differential cross 
section of e + e~ 77/) versus the invariant mass of a photon and p is drawn in Fig. H]as an 
example for the same compactness reason. 

TABLE I: The cross section of e + e~ -> -/VV at y/s = 10.58 GeV. 



final states 


invariant mass 


of VV (GeV/c 2 ) 


a(fb) 


7 J/i> p 


3.8 


- 5.2 


0.29 


7 ■VV' 4> 


4.1 


- 5.2 


0.022 



7P0 1.7-5.2 0.66 



TABLE II: The cross section of e+e~ -> 7VV at = 3.097 GeV. 



final states 


invariant mass 


of VV (GeV/c 2 ) 


a(fb) 


7 P 


1.7 


- 3.0 


9.68 


7 u; 


1.8 


- 3.0 


0.80 



IV. DISCUSSION 

In this article, we calculate the cross sections of e + e~ annihilation into one photon plus 
J/iftp, J/^fp-, or pcf) at y/s = 10.58 GeV; and one photon plus ujcp, pep, or 7p, 70, or •yu 
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FIG. 3: Differential cross section of e+e~ — > 77*7* —> jiftp at y^i = 3.097 GeV. 
TABLE III: The cross section of e+e" -> 77F at = 3.097 GeV. 



final states 


invariant mass of VV (GeV/c 2 ) 


a (pb) 


7 7 P 


0.7-3.0 


49.30 


7 7 


1.0-3.0 


6.75 


7 7 UJ 


0.8-3.0 


3.96 
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FIG. 4: Differential cross section of e + e~ — > 777* — > 77/5 at y^i = 3.097 GeV. 



at yfs = 3.097 GeV. The numerical results are shown in Tables fl] [TTl and IHIj and several 
differential cross sections are displayed as examples in Figs. EJ [31 and HI 

One essential component of our calculation is about the processes e + e~ — > (71)7(771)7*, 



here we only calculate e + e 



1)7(2)7* and e H 



(2)7(1)7* at leading order. However, 



it should be pointed out that our method is flexible and extensible, i.e., within the frame- 
work of our computation, the annihilation is a purely QED process, any number of photons 
and higher orders can be achieved without difficulty For any extension, the whole method is 
standard, only a few collateral modifications need to be carried out within the computation 
frame. More details are given in the following paragraph. We should also mention that, up 
to now, we have neglected some possible corrections from such as hadronic contributions, 
higher order loops, weak interaction, product photons interference, and the mass distribu- 
tions of the resonances. However, these corrections are small compared with the leading 
order contribution and far beyond the present experimental precision. Similarly, the above 
improvement can be achieved through just intuitive extensions of our present calculation. 

Here we take some interesting deductions from our calculation such as e + e~ — > VV and 
e + e~ — > 77(7). In addition to illustrate the flexibility of our method, these deductions 
also conform the validity of our computation by the consistencies between them and other 
literatures. First, it is obvious that when we reduce the number of photons from three 
to two, we are actually calculating e + e~ annihilation into two virtual photons. Then it is 
easy to get the cross sections with two-vector-meson final states such as e + e~ — ► uJ/ip, 
e + e~ — > p<p, and e + e~ — > pJ/ip. It turns out that the results from our calculation are the 
same as those in Ref. [l3[ ■ Note that the whole calculation is standard, only a few specific 
considerations are taken into during this reduction, such as the 2 — > 2 phase space and 
package caused" factors 



some 



141 ] which are different from those in the process 2 



Second, we can fix all the final state virtual photons to be real, then actually we calculate 
the cross sections of e + e~ — ► 777. Similar to the above case, here we need only do a few 
peripheral modifications, that include fixing all the photon-masses to zero since the photons 
are on-shell now, setting the number of polarization directions to be two instead of previous 
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three, and multiplying a factor 1/3! = 1/6 because now the final state contains just three 
identical bosons etc. Eventually our result of e + e~ annihilation into three real photons is 
consistent with both theoretical prediction [l5[ and experimental result [23[. Finally, with 
a further reduction when we only calculate e + e~ annihilation into two real photons, the 
analytic formula returns to the familiar (1 + cos 2 9)/ sin 2 9 form. All the reductions we 
discussed above are easily realized in our program and their consistencies with other studies 
should be considered as verifications of our method. 

Finally we want to mention two important features of our results. One is that if the phase 
space allowed, i.e. the CM energy is high enough for a specific final state, all the differential 
cross sections would show similar shapes. As what displayed in Figs. [2] and HI with the 
invariant mass varying from low to high, there would be a bump near the threshold followed 
by a flat part and then ended with a fast increase when the invariant mass approaches 
the CM energy, which is resulted from a very soft radiated photon. The other feature 
is, as we expected, these three-photon processes are 0(a) suppressed compared with the 
corresponding ISR processes where the hadron system is produced from a single photon. 
That means the backgrounds which arising directly from e + e~ — > jVV and e + e~ — > 'yyV are 
not essential at £?-factories in the study of the ISR processes, nor at BES in J / ip decays at 
current available statistics. However, accompanying with the upcoming super- factory and 
BESIII, these modes will be important in the near future when more accumulated luminosity 
is achieved and a better precision is expected in various analyses. 
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